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Abstract. Light sterile neutrinos might mix with the active ones and be copiously 
produced in the early Universe. In the present paper, a detailed multi-flavor analysis 
of sterile neutrino production is performed. Making some justified approximations 
allows us to consider not only neutrino interactions with the primeval medium and 
neutrino coherence breaking effects, but also oscillation effects arising from the presence 
of three light (mostly-active) neutrino states mixed with two heavier (mostly-sterile) 
states. First, we emphasize the underlying physics via an analytical description of 
sterile neutrino abundances that is valid for cases with small mixing between active 
and sterile neutrinos. Then, we study in detail the phenomenology of (3 + 2) sterile 
neutrino models in light of short-baseline oscillation data, including the LSND and 
MiniBooNE results. Finally, by using the information provided by this analysis, we 
obtain the expected sterile neutrino cosmological abundances and then contrast them 
with the most recent available data from Cosmic Microwave Background and Large 
Scale Structure observations. We conclude that (3 + 2) models are significantly more 
disfavored by the internal inconsistencies between sterile neutrino interpretations of 
appearance and disappearance short-baseline data themselves, rather than by the used 
cosmological data. 
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1. Introduction 

During the last several years the physics of neutrinos has achieved remarkable progress. 
The experiments with solar [H [21 El HI El [6], atmospheric [3 E], reactor [9], and also 
accelerator neutrinos [10j HH [12] , have provided compelling evidence for the existence 
of neutrino oscillations, implying non-zero neutrino masses. 

The present data require at least two large (6 I 12 and #23) angles in the neutrino 
mixing matrix and at least two mass squared differences, Am]; = m] — mf (where 
rn^'s are the neutrino masses), one driving the atmospheric (Am^) and the other one 
the solar (Am^) neutrino oscillations. The mixing angles 812 and 823 control the solar 
and the atmospheric neutrino oscillations, while the third angle #13 in the three neutrino 
mixing matrix is bound to be small by the data from the CHOOZ and Palo Verde reactor 
experiments [T3l [H] . Unfortunately, oscillation experiments only provide bounds on the 
neutrino mass squared differences, i.e. they are insensitive to the overall neutrino mass 
scale. 

Cosmology provides one of the means to tackle the absolute scale of neutrino 
masses. The Universe can be thus exploited as a new laboratory where to test 
neutrino masses and neutrino physics. The new accurate measurements of the Cosmic 
Microwave Background (CMB) anisotropy and polarization from satellite, balloon-borne 
and ground-based experiments have fully confirmed the predictions of the standard 
cosmological model (see e.g. [15J) and allow us to weigh neutrinos. Neutrinos can 
indeed play a relevant role in large scale structure formation and leave key signatures in 
several cosmological data sets. More specifically, the amount of primordial relativistic 
neutrinos change the epoch of the matter-radiation equality, leaving an imprint on both 
CMB anisotropies (through the so-called Integrated Sachs- Wolfe effect) and on structure 
formation, while non relativistic neutrinos in the recent Universe suppress the growth 
of matter density fluctuations and galaxy clustering (see e.g [16]). 

All these new observations have allowed us to place new bounds on neutrino physics 
from cosmology [TFJ HH HH [15l [20J, [22], [23] constraining the amount of energy density in 
neutrinos and their mass, in the framework of three active degenerate neutrinos, to be 
below ~ 0.25 eV at 95% C.L. However, the three neutrino scenario is a minimal scheme, 
inspired, following a theoretical prejudice, by a model which tries to resemble the three 
family structure of quarks and charged leptons. We now know neutrinos are different 
from the quark and charged lepton sectors and there is no fundamental symmetry in 
nature forcing a definite number of right-handed (sterile) neutrino species, as those are 
allowed in the Standard Model fermion content. 

Models with three active neutrinos plus an additional sterile one (3 + 1 models) were 
introduced [24J to explain simultaneously solar, atmospheric and LSND [25] data. The 
mass of a fourth, thermal, sterile neutrino has been constrained from cosmological data 
to be less than ~ 0.65 eV at 95% C.L. (see e.g. [19J). While cosmological constraints 
on the neutrino mass are certainly remarkable, one should be careful in taking those 
constraints blindly, since they are obtained in a model-dependent way. For instance, 
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the constraint on the mass of the fourth sterile neutrino presented in [19j is obtained 
under the assumption of a (3 + 1) scenario with a fully thermalized sterile neutrino. On 
the other hand, if different thermalization mechanisms are considered [26], |27[ [28], 29J, 
the cosmological constraints on neutrino masses could be relaxed or even drastically 
changed. 

More recently, scenarios with three active plus two sterile neutrinos (3 + 2 
models) [30] have been shown to provide a better fit to short baseline (SBL) data 
compared to (3 + 1) models [3TJ [32, [33], since they can, within a neutrino oscillation 
framework, bring into agreement the LSND v e excess in a beam [25] and the null 
results from other SBL oscillation experiments. The masses for the additional mostly- 
sterile neutrinos we are interested in are much larger than the light neutrino ones, 
m 4 ,m 5 3> mi,m 2 ,m 3 . Throughout this paper, we take the convention m 5 > m 4 . The 
flavor neutrino base v a , a = e, //, r, s,p, is related to the massive base z/j, i = 1,2, 3, 4, 5, 
through a 5 x 5 unitary matrix which we indicate as U: 

V a = XJuiVi . (1) 

Searches of sterile neutrinos indicate that the mixing between heavy neutrinos and active 
neutrinos and between light neutrinos and sterile ones is very small, so throughout our 
study we will take \U a ^, \U a ^\, \Uj S \, \Uj P \ <C 1, with a = e, \i, r and j = 1, 2, 3. 

The way in which these neutrinos shape the Universe depends crucially on the 
thermalization processes. The purpose of the present study is to test the cosmological 
viability of the (3 + 2) scenario in which sterile neutrinos are produced through 
oscillations |Ml 135] ESI [33 EHl [39l [40l HIl [4^ (see [43l HH [45] for earlier works), applying 
a multi-flavor treatment in the evolution of the system. We add in the analysis the most 
recent neutrino oscillation data, which includes the first MiniBooNE — > v e oscillation 
results. Finally, let us note that the presence of additional relativistic degrees of freedom 
at the epoch of Big Bang Nucleosynthesis (BBN) modifies significantly the production 
of helium, deuterium and lithium (a recent review is given in [46], see also [47]), so 
this observational data on the abundance of light elements allows one to put severe 
constraints on sterile neutrino parameters. However, in the present study we consider 
the contribution of sterile neutrinos to the energy density of the Universe and its effects 
on the CMB and on the matter power spectrum, but we do not include the data from 
BBN observations. 

The structure of the paper is as follows. In Sec. [2] we describe the sterile neutrino 
production. Although we do not attempt to solve the exact quantum kinetic equations 
for a (3+2) system, we have numerically computed the momentum-averaged counterpart 
with an approximation for the coherence-breaking terms. This amounts to solve an 
equation for the 5x5 (anti)neutrino density matrix in order to use these solutions in our 
Markov chain Monte Carlo (MCMC) analyses. We present as well an analytical approach 
which, under a number of approximations, describes the neutrino thermalization in a 
multi-family scenario, deriving very simple formulae that illustrate the evolution process 
in a very clear way. Sec. [3] describes the SBL oscillation data and the cosmological data 
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sets used in the MCMC analyses. We present our results in Sec. HJ discussing first the 
compatibility of oscillation and cosmological data and presenting afterwards a combined 
analysis of both data sets. We conclude in Sec. El 

2. Sterile Neutrino Production in the Multi-Flavor Approach 

As a case study of sterile neutrino cosmological abundance in multi-flavor schemes, 
we consider mixings among three active plus two sterile neutrino states, i.e., the so- 
called (3 + 2) sterile neutrino models. In the Early Universe, and in particular for 
temperatures above neutrino decoupling (T > 1 MeV), neutrinos are part of a gas of 
interacting multispecies of elementary particles. In this environment, sterile neutrinos 
can be produced via the oscillations of active neutrinos in presence of interactions with 
the thermal plasma (for foundational work see pi SHI H3 ES EH E2 EHl El EHl EH]). Due 
to the high temperatures, the neutrino mean free path is quite small and it is essential 
to take into account the breaking of coherence, using the density matrix formalism. The 
quantum operator for the density matrix of neutrinos is defined as p* = z/z/J\ The C- 
valued density matrix, p(E, T), is obtained by taking the average of the matrix element 
over the medium. Although the exact form of the terms which take into account the loss 
of coherence is quite complicated [51j E21 EH ESI EI], a good approximate treatment is 
obtained if the kinetic equations for p are written as in [58] (where we will omit writing 
the E and T dependence in what follows): 



Here, H m = UHqU^ is the free neutrino Hamiltonian in the flavor basis, which is 
obtained rotating the Hamiltonian in the mass basis Hq = diag(i?i, E2, E%, E4, E§) by 
the unitary mixing matrix U given by Eq. ([1]). For the temperatures under consideration, 
neutrinos are highly relativistic and the approximation E = p + m 2 /(2p), with p(~ E) 
the neutrino momentum, holds and will be applied. The effective potential V e g describes 
the interactions of neutrinos with the medium and is diagonal in the flavor basis. For 
negligible lepton asymmetry, its elements are approximately given by: 



where Gf is the Fermi coupling constant, T is the plasma temperature, E is the neutrino 
energy, a = 1/137 the fine structure constant. The constants C a are given by C e ~ 0.61 
and C^r ~ 0.17 (for T < m M ) [59 J, which are exact in the limit of thermal equilibrium 
of all leptons involved, and sufficiently accurate for our purposes. The coherence loss 
in the evolution of the non-diagonal terms is given by p = —Tp [131 EH 001 El] and a 
simplified but accurate way to describe the neutrino production and destruction is given 
by the anticommutator term in Eq. (T5]) [58J. The damping factor T is diagonal in the 
flavor basis, T = diag(r e , T M , T T , 0, 0), and its elements, taking the total scattering rate, 
in the Boltzmann approximation |4Q| . are given by: 



p = i [H m + Kff, p] - {r, (p - p cq )} • 



(2) 



C a G 2 F T 4 E/a , 



(3) 



90C(3) 
7vr 4 




(4) 
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with a = e, //, r. The coefficients g a have been computed numerically and are given by 
g e ~ 3.6, = g T ~ 2.5 (60]. The diagonal matrix p eq = diag(p eg ) = I(exp(E/T) + 
with / the identity matrix, is the equilibrium value of the density matrix. 

For the purpose of this work we can assume that, in the range of temperatures 
under interest, the time dependence of the temperature scales as T ~ —HT, where H 
is the Hubble expansion rate. Using the fact that 

(5) 

E_ 

T 



we can rewrite Eq. 

/ dp 




dT I E HT 



[i[H m + V cS ,p}-{T,(p-p cq )}) , (6) 



where 



H(T) = J^IL . (7) 

v ' V 45 mpi y ' 

Here, mpi is the Planck mass, while g+ is the number or relativistic degrees of freedom. 
In the range of temperatures where all the processes we are interested in take place, 
g* changes very little. Hence, we have neglected the temperature dependence of g* in 
the time dependence of the temperature, although we have taken it into account in the 
expression of the Hubble rate. 

The MCMC analyses presented in Sec. |3] use the solutions to the approximate 
evolution equation, Eq. ([6]), for the 5x5 (anti) neutrino density matrix p. Given the 
hermiticity of the C-valued density matrix, this amounts to solving a set of 30 coupled 
first-order differential equations. Thus, in order to reduce the computing load, to solve 
Eq. ([6]), we only consider monoenergetic neutrinos at any given temperature, taking an 
average value of E = 7£(4)T/(2£(3)) ~ 3.15 T. In general, if the oscillations are faster 
than the Universe expansion rate as is the case, this is a justified approximation [56J. 
On the other hand, we take thermal abundances for the mostly-active neutrino triplet 
and zero abundances for the mostly-sterile states, as initial conditions at a temperature 
of T ~ 100 MeV, and follow the temperature evolution until T = 1 MeV0. As a further 
approximation, we assume active neutrinos keep perfect thermal distribution throughout 
the evolution. This is justified as any active neutrino population depleted due to 
oscillations will be immediately repopulated thanks to its interaction with the plasma. 
Furthermore, if there is no neutrino-antineutrino asymmetry as we are assuming, active- 
active oscillations must have no effect at all HB1 [561. 



In the following, Sec. 12.11 we present an analytical description of the solutions to 
Eq. ([6]), providing formulas for sterile neutrino production in the Early Universe that 
are very good approximations for sterile neutrino scenarios with small mixings among 
active and heavy neutrinos (and among sterile and light neutrinos). We will use these 
results to show the effects of different approximations assumed when solving numerically 
the problem. 

| Note that post weak decoupling effects can only be important when there is a large lepton number [61j . 
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2.1. Analytical Description 

The production of sterile neutrinos can be computed by solving Eq. (j2j) and in particular 
by following the evolution of the diagonal sterile terms p ss and p pp . In the present 
analytical derivation, we are interested in values of the heavy sterile masses such that 
Am^j, Amy ^> Am^, Am^ and small mixing between the active and heavy sectors and 
between the sterile and light sectors. In this case, sterile neutrinos never reach complete 
thermalization and their abundances are much lower than the equilibrium one. 

As mentioned above, the resolution of Eq. ([2} implies solving a system of 30 coupled 
first-order differential equations. However as the typical frequency for oscillations 
between active (u a , a = e, [A, t) and sterile (uh, h — s,p) neutrinos, given by 



U a h 



is much larger than the expansion rate of the Universe H in the epoch when sterile 
neutrino production takes place, we can use the static approximation [62"1 [57]: 

Pah = , (9) 

Ps P = o . (10) 

All the off-diagonal terms can be found by solving this system of 20 coupled linear 
equations. Hence, this approximation reduces the problem to the resolution of a system 
of 5 coupled first-order differential equations, involving only the diagonal elements of 
the density matrix, which is given by 



Pact 2 

where we have defined: 



J2 Im ( H apPap) ~ Fa (paa ~ Pcq) 



a,p = e,p,r,s,p (11) 



(12) 

i 

with Kfr,a given by Eq. El and V eS>s , V eSjP = 0. 

For the small mixing case considered, production of sterile neutrinos never reaches 
thermal abundances, i.e., p ss ,P PP "C p eq , and the effects on the active neutrino 
distribution are very small and therefore we can use the approximation p aa ~ p cq [§|. 
This simplifies further the problem and, using the fact that H a p H sp , H ss , H pp holds 
in this case, with (3 = e,p,r,s,p, the remaining three systems of four linear equations 
corresponding to Eqs. (El), one system for each active neutrino, can be solved to find: 
^ H as (H aa — H pp + zr o ) + H ap H* p 

Pas — Peg 777 7J w ~ r j ,„ , [i-O) 



(H aa — H ss + iT a ) (H aa — Hpp + zr a ) — \H. 



sp I 



H ap (H a a H ss -\- ir a ) -|- H a sH sp , 

Pap — PeqJT} 7} -t \ / tj tj , 7F \ I tt 12 ' v > 



(H aa — H ss + iT a ) (H aa — Hpp + iT a ) — \H. 



sp I 



§ Note that we will also use this approximation for larger mixings, as those studied below. Even in 
that case, the distortion of the active neutrino distribution is small enough for our purposes and is not 
expected to alter substantially our results. 
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Although p sp turns out to be subdominant with respect to p aS} it plays a relevant role 
in the evolution of p ss and p pp and cannot be neglected. From Eq. fflQj) . we get: 

Ea (PlsHap - PapH* s + H sp (p ss - p pp ) 



Psp — Peq' 



H ss — H 



(15) 



pp 



In the approximations we made, the evolution of p ss and Ppp CcLll be explicitly given 



by: 



Ps 



Ppp 



-2 Im (j2 PasKs - PspH* p ^j , 

-2 Im (^2 PapH* ap + Ps P H*p^j . 



(16) 



;i7) 



By substituting the expressions for p ah in Eqs. ( TTBl and ffT4l) and p sp in Eq. ffl5|) . one 
can obtain the two differential equations for p ss and p pp which can then be solved in 
order to find the sterile neutrino abundance in the Early Universe. It is interesting to 
notice that the terms proportional to p ss and p pp in the imaginary part in the r.h.s. of 
Eqs. (1161) and ffTTl) are real as they are proportional to \H sp \ 2 and therefore vanish. This 
implies that the system is now given by two uncoupled differential equations which can 
be solved either numerical or analytically with some approximations. 

For a real mixing matrix, which will hold in the numerical study we will perform 
later, the results are simplified and the evolution of p ss and p pp depends only on the 
imaginary part of the off-diagonal elements of the density matrix p^p, h = s,p, as: 



Phh 



-2 



J2 H h0 Im(p 0h ) 

which we will use throughout the remaining of this section. For compactness of notation, 
we define 



I a p = lm(p a p) = -lm(pp a ) . 
From Eqs. ( f!3l) and f Cbfj) . we can find the explicit form of I as and I ap '- 



(19) 



I ~ 

± as — 



{H aa - Hppf + T 2 a + H ; 



D 



+ 



[{H aa — H ss ) + {H aa — Hpp)} H sp H, 



ap 



I ~ 

± ap — 



[H n 



D 

H ss f +TI + H 



(20) 



H, 



ap 



D 



+ 



[{Hga — Hpp) + {Hgg — H SS )] H Sp Hg 

D 



Peq 



where 



D 



(Hgg — Hpp) (Hgg ~ Hg 



H 



up 



+ 



{Hgg — H ss ) 2 + (Hgg — Hpp) 2 + 2H 



sp 



+ r 



(21) 

(22) 
(23) 
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Finally, from Eq. (IT5|) . the element I sp can be further expressed in terms of the six I a ^, 

Isp — ~Tj 77 ^ ] {H a pI(x S ~\~ H as I a p) . (24) 

Note that, as already discussed for the general complex case, I sp turns out to be 
subdominant with respect to I a h, as it is suppressed by H ah /(H SS — H pp ). However, 
as it is clear from Eq. ([TBI) , it can play a relevant role in the evolution of p ss and p pp . 
In that equation, whereas the prefactor of I a h is H a h, that of I sp is H sp , which can 
compensate the smallness of I sp with respect to I a h- Hence, whenever there is mixing 
between the two sterile states, I sp must be taken into account. It is important to note 
that I as , I as and I as are proportional to T a , indicating that the production of sterile 
neutrinos in the Early Universe is due to the breaking of coherence in the evolution of 
active neutrinos. 

Using the approximation that the time-dependence of the temperature scales as 
f ~ -HT and Eq. ©, we can rewrite Eqs. (lTg]) - pE]) as: 



HT f^l -2Y J H as I as -2H sp I sp , (25) 




HT ~2Y,H ap I ap + 2H sp I sp . (26) 



This is the system of two uncoupled first-order differential equations which has to be 
solved. Indeed obtaining the final abundance of sterile neutrinos can be easily performed 
numerically. 

In addition, in the case of no mixing in the sterile sector as we are considering in 
this work, i.e., U s5 , U P 4 <C 1 (H sp ~ 0), Eqs. (|20l) and (|2~T!) further simplify (with h = s 
or p) to: 

T _ H ah 

(H aa — H hh ) + r 

and phh is given from Eqs. ( 1251) and ( 1261) by 



ah — f TT TT N 2 , t^2 ^ a P e 1 ' (27) 



n roo 1 



Peq J Tdcc HT ^ (H aa - H hh Y + VI 

where E/T is kept constant. Here, T dec is the decoupling temperature of the active 
neutrinos, which is of few MeV, depending on flavor. The term T 2 a in the denominator 
can be neglected as this gives a small relative correction to the final value of the sterile 
neutrino distribution of order (T a /V a ) 2 ~ 10~ 4 . Using this, we can rewrite Eq. (128]) . 
keeping only the dominant terms and substituting numerical factors, as 

0.053, f^y % EfcDj r 773 * _,(») 



Peq V eV ) y/g+ a J Jx d 



Am ii 



^r^- + 3.7x 10- 8 C a a; V 



where x = (T/MeV) 3 (a; dcc = (T dec /MeV) 3 ), y = E/T and j = 4 (5) for h = s (p). This 
is the same integral which appears when solving the simplified two-neutrino case 
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A further simplification can be applied by taking into account that, for the values 
of Am| x and Am^ we are interested in 0{1 eV 2 )), the maximum of production of 

1 /6 

sterile neutrinos happens at T max ~ 13 MeV (Am^/eV 2 ) [H [36j EH H2J [58] and thus 
^max ^ ^dec) so the lower limit of integration can be safely put to 0. The relative error 
introduced is of order T^GpE /TyJCJ (a Amf x ) ~ lCT 3 ^! eVVAm 2 !. With these 
approximations the integral can be performed analytically and the final sterile neutrino 
distribution is given by: 

^^6.6x 10- 3 ,S E -% (H) 2 • (30) 
p eq V eV 2 V v 7 ^ UO- 2 / 

where we have taken = 10.75 and £7jy — 1. It should be noticed that for large values 

of U a j (~ 10~ 2 ), the sterile neutrino distribution approaches the equilibrium value and 

our results are not valid anymore. 

As Phh{E,T) has the same functional form as p eq (E,T), it is straightforward to 

obtain the contribution to the energy density of the Universe of each of the heavy 

states, which is simply given by 

We reiterate that the results in Eqs. (130!) and (I3TI) are only valid for p ss , p pp -C p eq . It 
is straightforward to recover the well-known two-neutrino results from these expressions, 
which as one could expect, indicate that in the case of no mixing between the sterile 
states, the calculation of each of the heavy neutrino abundances reduces to adding the 
individual contributions due to their mixing with each of the active neutrinos. 

2.2. Accuracy of approximations 

As mentioned above, we do not attempt to solve the exact quantum kinetic equations for 
a (3 + 2) system, but instead numerically find the momentum-averaged solutions for the 
simplified counterpart given by Eq. This approach describes the loss of coherence in 
an approximate way which, along with the assumption of monochromaticity, allows us to 
reduce the integro-differential system into a system of first-order differential equations. 

It is well known that the averaging over momentum is a justified approximation 
if the relative oscillation phases can be ignored, as is the case if oscillations are not 
slower than the expansion rate of the Universe [56]. On the other hand, it has also been 
shown [58] that a reasonably accurate description of the coherence loss can be achieved 
by making use of the anticommutator in Eq. ([6]). 

Keeping in mind that our starting point is not the system of exact equations but an 
approximate one, we discuss some further approximations assumed in order to be able 
to fully solve the system of 30 coupled first-order differential equations and study what 
is their effect on the final result. In order to do so, we consider the analytic calculation 
described above for the small mixing case. 

Among the different approximations used, there is the monochromaticity of the 
neutrinos, assuming that active neutrinos are always in equilibrium, neglecting the 




T [MeV] 

Figure 1. Ratios p ss /Peq and p pp /p e q as a function of the temperature as given by 
Eq. (|28p with all approximations (see text) and when only some of these approximations 
are implemented. The black (red) solid lines refer to the case when all approximations 
are taken to obtain p ss (p pp ). The blue (orange) lines represent solutions to p ss (p pp ) 
with some approximations taken more accurately; dashed lines refer to the case when 
the momentum spread is taken into account, dot-dashed lines to the case when the 
exact form of the effective potential is considered and solid lines represent the case when 
both the momentum spread and the exact potential are assumed. The parameters 
we have taken are: U e 4 = = 2 x 1CT 3 , U e § = U P 5 = 10 -3 , m± = 1 eV and 
m 5 = y/W eV. 

temperature dependence of the number of degrees of freedom and using an 
approximate expression for the effective potential. 

We find that the most important effect when sterile neutrinos are close to 
equilibrium is due to neglecting the change of the number of degrees of freedom as 
they get produced. Obviously, this approximation has a negligible effect in the small 
mixing case. Surprisingly, however, the effect of not taking into account the departure 
from equilibrium of the leptons as the temperature decreases, i.e., using an approximate 
expression for the effective potential, is comparable to the previous one even in the case 
of small mixing. On the other hand, we have explicitly checked for various cases that, 
as expected, solving for an averaged momentum introduces a very small error. 

To conclude, these approximations tend to overestimate the final abundance, 
although none of them induces an error larger than ~ 10%. This is depicted in Fig. [TJ 
where we show the effect of two of these approximations in the case of small mixings by 
solving Eq. (1281) . The black (red) solid lines represent the solution for p ss (p pp ) when 
all approximations are implemented. By removing these approximations, we show their 
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Figure 2. Flavor composition of neutrino mass eigenstates in (3 + 2) sterile neutrino 
models. The hatched rectangles indicate active flavor content (electron, muon, tau, 
respectively, from left to right), and the empty rectangles indicate sterile flavor content. 

effect on the final solution. The blue (orange) lines are the solutions to p ss (p pp ) when 
some of the approximations are not taken. The dashed lines refer to the case when the 
momentum dependence is taken into account, the dot-dashed lines represent the case 
when the exact form of the effective potential is considered and the solid lines depict 
the case when both the momentum spread and the exact form of the potential are 
assumed. As one can see the effects are relatively small and hence we are confident that 
the results we obtain are accurate enough for our purposes. Finally, let us note that 
we have also checked the perfect agreement (in this limit) between the semi-analytic 
solutions of Eq. (1281) and the numerical ones of Eq. ([6]). 

3. Data Samples and Analysis Procedure 

We now consider the signatures of the heavy mostly-sterile neutrinos in neutrino 
oscillation experiments and we study the sterile neutrino production in the Early 
Universe with a detailed numerical simulation. The results of this simulation are valid 
both for the p ss , p pp <C p eq case discussed above, as well as for the case where sterile 
neutrinos reach thermal abundances (p ss ,p pp ~ p eq ). 

A schematic representation of (3 + 2) models considered in this analysis is shown in 
Fig. 2. We take a simplified form for the 5x5 mixing matrix, requiring no mixing with 
the r-neutrino sector. For the two mostly-sterile states, we allow for variable electron and 
muon flavor content ( | f/ e4 1 2 , | ^7^4 1 2 , |cV e5 | 2 , |t/ M5 | 2 ), and variable neutrino masses m 4 and 
ms, to be constrained by our analysis described below. Concerning the three mostly- 
active mass eigenstates, we fix the mass and mixing parameters within the regions 
currently allowed by solar, reactor, atmospheric, and accelerator long-baseline neutrino 
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experiments. Specifically, we take m\ = 0, m 2 = y/8 x 10 -5 eV, m 3 = x 10~ 3 eV 
for the masses [63], and the upper left 3x3 block form given in the following 5x5 mixing 
matrix: 



/ 



0.81 
-0.51 
0.28 







0.55 
0.51 
-0.67 








0.70 
0.70 







±\u e4 \ 
±\uva\ 



1 





±\u e5 \ 




1 



(32) 



where the index a = e,fi,r,s,p runs over flavor states, the index i = 1,...,5 over 
mass states, and we indicate the two sterile states with the symbols s, p, respectively. 
The following simplifying assumptions were therefore made: the lightest neutrino mass 
eigenstate is assumed to be massless; the normal hierarchy for the mostly-active triplet 
is taken; the mixing angle $13 is taken to be zero, and $12, $23 are fixed to their current 
best-fit values [63]; all six Dirac CP- violating phases allowed for 5 neutrino species are 
assumed to be zero; zero r content in the 4th and 5th mass eigenstates is assumed; 
and no mixing between the sterile states (U s c, = U P 4 = 0). All but maybe the first 
assumption are expected to have a negligible effect; relaxing the mi condition would 
impose stricter bounds on (3 + 2) models with respect to what reported in the following. 
We also note that, while do not know explicitly enforce exact unitarity conditions in 
Eq. (1321) [{J, this is approximately ensured by the smallness of the (|?7 e4 | 2 , |C^4| 2 , |^es| 2 , 
I 1 2 ) elements, providing a sufficiently accurate description. In summary, our analysis 
of SBL and cosmological data allows us to constrain six parameters of (3 + 2) models: 



the masses 777,4, and the mixing matrix elements |t/ e 4| 2 , \U^\ 
following, combinations of these six parameters will be shown. 



\U. 



e5 I 



\u, 



/j,5 1 



In the 



3.1. Short- Baseline Oscillation Data 

Our analysis of SBL data follows closely what done in J3TJ [32]. We consider 
neutrino and antineutrino data samples for three different oscillations channels: 
muon-to-electron neutrino transitions, electron neutrino disappearance, and muon 
neutrino disappearance. Muon-to-electron neutrino data sets considered come from 
the LSND [25J, KARMEN2 [M], and NOMAD [65] experiments. With respect to 
previous analyses, we also include recent data [66J on the z/ M — > v e search by the 
MiniBooNE experiment [67] . Constraints on electron neutrino disappearance are based 
on data from the Bugey [68] and CHOOZ [13] reactor experiments, while data from 
the accelerator-based CDHS |69j and CCFR [70] experiments were used for muon 
neutrino disappearance. In addition, muon neutrino disappearance is constrained also 
by atmospheric (7J[8] and long-baseline [10] data, as done in [71]. The SBL experimental 
input used is summarized in Tab. [Q 

Given the baselines and neutrino energies available to the short-baseline 
experiments described above, it is possible to constrain oscillation frequencies in the 

|| Note that not even the 3x3 active neutrino block is taken as exactly unitary. 
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MiniBooNE 






0.2-1.1 


8 


NOMAD 






0.002-0.3 


30 


Bugey 


v e -> 




2-50 


60 


CHOOZ 






100-400 


14 


CCFR 


v„ -> 




0.004-0.03 


18 


CDHS 






0.02-1.5 


15 



Table 1. Summary of SBL data used in this analysis. 



Data Set 


m 4 (eV) 


m 5 (eV) 


u ei 




u e5 


u, 5 


SBL-only 


0.96 


5.0 


0.12 


0.15 


0.59 x 10" 1 


0.16 


SBL+cosmo 


0.68 


0.95 


-0.37 x lO" 1 


0.77 x 10~ 2 


0.13 


0.19 



Table 2. Summary of (3 + 2) model best-fit mass and mixing parameters. Top row 
using SBL data only and bottom row for a combined fit of both SBL and cosmological 
data. 



LSND allowed region, i.e., Am 2 ~ 1 eV 2 , as well as the corresponding oscillation 
amplitudes involving electron and muon neutrino flavors. Data are fitted to the (3 + 2) 
appearance and disappearance oscillation probabilities applicable for short baseline 
experiments and described in [3H [32], allowing to extract allowed regions in the six- 
dimensional mass and mixing parameter space considered in this analysis. Parameter 
values that best describe SBL data are given in the top row of Tab. [2j 

3.2. Cosmological Data and Model 

With the advent of the modern CMB datasets, such as those coming from experiments 
like WMAP [72] , Boomerang [73] and ACBAR [74] data, combined with modern galaxy 
surveys, for example SDSS and 2dF surveys, cosmology has entered a precision era in 
which not only basic cosmological parameters can be constrained, but also fundamental 
physical parameters that enter the theory. 

Cosmological neutrinos have a profound impact on cosmology since they change 
the expansion history of the Universe and affect the growth of perturbations [75J. 
Cosmological probes are sensitive to the number of neutrinos, their relative masses and 
abundances. In order to get theoretical predictions for cosmological datasets, Boltzmann 
codes must be used that follow the cosmological perturbation evolution in the linear 
regime, accurately accounting, in addition to the standard cosmological constituents, 
for each neutrino species. In this work we use the popular package cosmomc |76j . which 
has been adapted to work with more than one massive neutrino species. 

Our basic cosmological model is the minimal inflationary ACDM model that is 
consistent with most standard cosmological probes. In addition to neutrino parameters 
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discussed below, its parameters are the baryon density expressed as a fraction of the 
critical density multiplied by h 2 , = Vl^ti 1 , where h = ifo/100km/s/Mpc is the reduced 
Hubble constant H , same for the cold dark matter Udm = ^dmh 2 , the ratio of sound 
horizon to the angular diameter distance 9, which is essentially a proxy for the Hubble 
constant, the optical depth to the last scattering r, logarithm of amplitude log A and 
spectral index of primordial fluctuations n s . These basic parameters are enough to 
provide a good fit to all available data and set other parameters. For example, the dark 
energy density in a flat Universe is simply given by = 1 — — ^W- 

The cosmomc code has been expanded to include 5 species of neutrinos. We assume 
that their masses are given by 

mi = (0, 9 x l(T 3 eV, 5 x l(T 2 eV, m 4 , m 5 ) (33) 

and abundances by 

fli = (l,l,l,o 4 ,o 5 ). (34) 

This results in four extra parameters, m 4i5 and a 4i5 , on top of those that set the 
basic cosmological model. In principle, indices 4 and 5 are interchangeable and therefore 
likelihood space is symmetric with respect to this change. We do not attempt to 
change the parametrization to account for that, but instead use the symmetry as an 
additional check for chain convergence. We use a flat prior between zero and unity for 
the abundance parameter Oj. Since at zero abundance, the mass can take any value, we 
limit the masses to be less than 20 eV. In practice, this upper limit is not relevant for 
the abundances and masses of interest. 

On the data side, we start with a conservative compendium of cosmological datasets. 
These include the standard cosmic microwave background data, namely WMAP 5- 
year data [15], [72], Boomerang 03 data [73], the latest ACBAR data [H] and the 
VSA data [77] . In addition we use the data on the matter power spectrum from the 
spectroscopic survey of Luminous Red Galaxies (LRGs) from the SDSS survey [78J. To 
constrain the basic model, we also use the constraints coming from the latest compilation 
of supernovae [79]. Finally, we also apply the prior on the reduced Hubble constant of 
h = 0.72 ± 0.08 from the Hubble key project [HQ]. 

The basic cosmological model and the basic parameter set is what we use throughout 
the paper and it represents conservative constraints from the minimal cosmological 
model. The check the sensitivity of our constraints on the model used, we also calculate 
constraints on an expanded model, in which we also vary the total energy density of 
the Universe, or equivalently the curvature, parametrized by fi^ = 1 — J2&>i, where the 
index i runs over the various components of the Universe and the parameter w (pressure 
over density) that describes the equation of state of the dark energy. In the concordance 
model, these two parameters have values zero and —1, respectively. On the other hand, 
we also use cosmological data more aggressively and use the Lyman-a forest data from 
the SDSS quasar sample (see e.g. [20]). This dataset provides very good constraints 
on the neutrino properties and has no known systematics at the moment, but it is an 
extremely difficult experiment that has proven to be somewhat controversial. 
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Figure 3. Evolution of sterile neutrino abundances, normalized to thermal equilibrium 
abundances peq, as a function of temperature T. The dashed red (solid blue) line 
indicate p ss /peq (Ppp/peq)- The case shown corresponds to the (3 + 2) model that 
describes SBL-only data best. 

3.3. Analysis Procedure 

In practice, we proceed as follows. First, we generate via a MCMC method a large 
number of (3 + 2) models that give a potentially viable description of SBL data. 
Specifically, in this analysis we use about 1.6 x 10 4 models characterized by a SBL 
goodness-of-fit within about 20 x 2 units of the SBL best-fit. Second, for each one of 
those models, we determine the cosmological sterile neutrino abundances as described 
in Sec. |2j As one example, we show in Fig. [3] the behavior of the sterile neutrino 
abundances, normalized to the abundance at thermal equilibrium, as obtained by our 
analysis for the (3 + 2) model that describes SBL data best (see top row of Tab. [2]). 
Fig. [3] shows the evolution of p ss /p cq and p pp /p eq as the temperature of the thermal 
bath decreases, starting with null abundances at T = 40 MeV (right side of plot). For 
the SBL best-fit model, one can see from Fig. [3] that both sterile neutrino states reach 
thermal abundances before decoupling. In other words, in this case, their masses fully 
contribute to the sterile neutrino matter density (Q s + Q p )h 2 . Third, the viability for 
each of the models considered is further assessed by contrasting the predicted sterile 
neutrino abundances with cosmological observations, as explained next. 

In order to combine the short-baseline oscillation data with the cosmological 
constraints, we first run extremely deep MCMC with about 70 - 100 thousands accepted 
samples. Each sample is sampled from the combined probability density in the N- 
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Figure 4. A% 2 profiles as a function of the sterile neutrino matter density (il s +fl p )h 2 . 
The red (blue) line indicates the case where cosmological-only (SBL-only) data are 
fitted. Three scenarios are shown for cosmological data fits: default data set and 
model (solid), data set including Lyman-a forest data and default model (dashed), 
default data set and cosmological model with free w and f2fc (dotted). The horizontal 
dotted lines define the 90% and 99% confidence level regions (1 dof). 



dimensional parameter space. The density of samples, projected to the 4-dimensional 
hyper-cube of parameters (7714,7715,04,05), thus corresponds to the probability density 
for those parameters, marginalized over the remaining cosmological parameters, which 
are, for the sake of this work, nuisance parameters. The relative probability of any model 
with some finite 4th and 5th neutrino masses and abundances relative to the standard 
minimal 3-neutrino model is calculated by comparing the relative densities of samples 
at the considered model position in the sample to the densities in the (04, a 5 ) = (0, 0) 
corner of the parameter space. We calculate the local densities of samples in spheres 
(which might be cut by the edge of the parameter space, but this is taken into account). 
These are then converted to the likelihood ratios, which are in turn converted to the 
effective A% 2 . This procedure allows us to importance-sample the MCMC from the 
short-baseline oscillation data in Sec. 13.11 

4. Results 

Two types of results are presented. First, in Sec. 14. 1[ we quantitatively address the 
compatibility of the SBL (including LSND and MiniBooNE) and cosmological data sets. 
Second, in Sec. 14.21 we perform a combined analysis of SBL and cosmological data. In 
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Figure 5. Allowed regions in (Am^d^ 2 + P^] 2 ), Aml^t/eg] 2 + |£/p 5 | 2 )) space, 
for SBL-only data (filled blue regions) and cosmology-only data (red contours). Light 
(dark) colors correspond to 90% (99%) confidence level regions (2 dof). The yellow 
star indicates the SBL-only best-fit point. 



both cases, the analysis is carried out under a (3 + 2) sterile neutrino hypothesis (see 
Sec. |3]) and using the multi-flavor description of cosmological sterile neutrino abundances 
(see Sec. H]). 



4-1. Data Sets Compatibility 

Fig. H] shows how the \ 2 f° r the SBL and cosmology data sets separately vary as a 
function of the sterile neutrino matter density (Q s + Q p )h 2 . The x 2 profiles for both 
data sets are shown relative to their respective best-fit \ 2 values. Cosmological data 
prefers a small (if non-zero) value for (Q s + Q p )h 2 , while SBL data is not consistent with 
a null (Q s + Q p )h 2 value, since the LSND v e excess cannot be interpreted in terms of 
neutrino oscillations in this case. Quantitatively, we find the following allowed intervals, 
with confidence levels for 1 degree of freedom (dof) given: 

• Cosmology: (Q 8 + tt p )h 2 <0.018 (90% CL), (Q s + Vt p )h 2 <0.028 (99% CL) , 

• SBL: 0.055< (fi, + Vt p )h 2 <0.081 (90% CL), 0.013< (fl, + Vt p )h 2 <0.097 (99% CL) 

In other words, no overlap is found at 90% CL in the (Q s + Q p )h 2 ranges allowed by the 
two data sets, while overlap exists at 99% CL. 

The previous numbers refer to our default cosmological data set and default 
cosmological model described in Sec. 13.21 the one used as reference for all following 
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Figure 6. Evolution of sterile neutrino abundances, normalized to thermal equilibrium 
abundances peq, as a function of temperature T. The dashed red (solid blue) line 
indicate p S s/peq (Ppp/peq)- The case shown corresponds to the (3 + 2) model that 
describes all data (SBL plus cosmology) best. 

results and figures. In Fig. HJ this case corresponds to the solid red line. We have also 
investigated the impact of varying the cosmological data set and cosmological model. 
The dashed red line in Fig. 0] uses an extended data set including Lyman-a forest data, 
providing a tighter cosmological constraint. On the other hand, the dotted red line 
in Fig. H] quantifies by how much the cosmological constraint is relaxed by assuming a 
cosmological model with free w and (see Sec. I3.2I) . 

Fig. O conveys similar compatibility information as the one shown in Fig. HJ but 
expressed in terms of particle physics rather than cosmological parameters. Two- 
dimensional allowed regions in (Am^dcT^ 2 + (f/^j 2 ), Am^i(\U e 5\ 2 + \U^\ 2 )) are shown, 
separately for the SBL and cosmological data sets. As the neutrino masses and/or 
the electron plus muon flavor of the two heavy (mostly-sterile) states increases, the 
heavy neutrino matter density (Q s + Q p )h 2 tends to increase, either because of the 
larger neutrino mass or because of a more fully thermalized abundance. Overlap is 
found between the regions allowed by SBL and cosmology only at 99%, but not at 90%, 
confidence level (2 dof). Given the cosmological constraints on (Q s + Q p )h 2 in Fig. HJ we 
note that the corresponding bounds on Am -j J2 a U 2 j, with j = 4, 5 and shown in Fig. [5J 
are in this case about two orders of magnitude less stringent compared to what would 
have been expected from Eq. [3TJ As previously mentioned, the reason is that Eq. [31] is 
only applicable for p ss ,p pp p eq and is not generally valid in this case. 
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Figure 7. Allowed regions in (Am^d^l 2 + P^] 2 ), Am^d^ 2 + \U^\ 2 )) space, 
for all data (SBL plus cosmology). Light (dark) blue regions correspond to 90% (99%) 
confidence level (2 dof). The yellow star indicates the global best-fit point. 



As an additional statistical test for compatibility, we consider the "parameter 
goodness-of-fit" (PG) defined in [ST]. The test allows to quantify how likely it is that 
SBL and cosmological results arise from the same underlying (3 + 2) sterile neutrino 
model. This method alleviates a problem affecting goodness-of-fit tests based on 
absolute x 2 values, namely that a possible disagreement between the two data sets 
is diluted by data points which are insensitive to the parameters that are common to 
both data sets. The number of parameters common to both data sets is 4 in our case, 
since cosmological data only depend on (777,4, 777,5, 04, a 5 ), where the sterile neutrino 
abundances 04 and a 5 defined in Sec. I3.2I are complicated functions of all mass and 
mixing parameters. This is because, in general, we cannot assume the sterile neutrinos 
to thermalize, and therefore the mixing matrix elements (in addition to the neutrino 
masses) play a role in the cosmological case also. The test is based on the statistic 



Xpg 



X PG , 



Xpg, sbl = (XsBiJaii min - (Xsbl)sbl min are the (positive) differences for the cosmology 
and SBL x 2 values obtained by minimizing the combined (cosmology plus SBL) x 2 
function discussed in Sec. I4.2I minus the x 2 values that best fit the two individual 
data sets. In our case, we obtain Xpg, cosmo = 1-71 and Xpg, sbl = 5.81, yielding 
Xpg = 7.52. The goodness of fit parameter corresponding to this Xpg value and four 
common parameters gives a 11.1% compatibility between the two data sets. 



+ Xpg, sbl, where Xpg, 



= (x 



cosmo /all min 
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Figure 8. Allowed regions in (Am| l7 Am^) space, for all data (SBL plus cosmology). 
Light (dark) blue regions correspond to 90% (99%) confidence level (2 dof). The dotted 
line indicates Am^ = Am^. The conventional choice Am^ > Am^ is adopted in 
this analysis. The yellow star indicates the global best-fit point. 



4-2. Joint Analysis 

Assuming that SBL and cosmological measurements arise from the same underlying 
(3 + 2) sterile neutrino model, we have also performed a combined analysis, where the 
function to be minimized is the sum of the SBL and cosmology \ 2 functions. The 
best-fit parameters for this combined analysis are given in the bottom row of Tab. [2j 
As expected, the inclusion of cosmological data tends to prefer smaller sterile neutrino 
masses, compared to the SBL-only case. In addition, the best-fit mixing matrix elements 
are such as to alter the thermalization of the lowest mass sterile state. This is illustrated 
in Fig. [6l which is the equivalent of Fig. [3] for the best-fit SBL+cosmological model, 
rather than for the SBL-only analysis. In this case, the fourth mass state decouples 
with a density of only about 44% of the thermal one. 

Allowed regions in (Am^d^l 2 + |^| 2 )' ^ m 5i(l^es| 2 + |^s| 2 )) space from the 
combined analysis are shown in Fig. [7J As expected, the allowed region in the top-right 
portion of Fig. [5] disappears with the inclusion of cosmological data in the fit, reducing 
the parameter space allowed to the overlap region between the SBL and cosmological 
regions. 

Figs. M and [9] show additional two-dimensional projections of the six- dimensional 
parameter space allowed by the combined analysis. Fig. [8] illustrates allowed values for 
the fourth and fifth mass states of (3 + 2) models, and Fig. [9] their total active flavor 
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Figure 9. Allowed regions in (|t/ e4 | 2 + W^ 2 , \U e 5\ 2 + \U l _ t5 \ 2 ) space, for all data (SBL 
plus cosmology). Light (dark) blue regions correspond to 90% (99%) confidence level 
(2 dof). The yellow star indicates the global best-fit point. 



content (electron plus muon). In these figures, we can see that favored (3 + 2) sterile 
neutrino models are characterized by one mostly-sterile state that is highly constrained, 
with Arrij! ~ 1 eV 2 , and electron plus muon flavor content of about 5%. On the 
other hand, mass and mixing parameters for the second mostly-sterile state are far less 
constrained, with mass splittings in the range 0.1-10 eV 2 , and electron plus muon flavor 
content in the range 10~ 5 -5 x 10 -2 . 

5. Conclusions 

Right-handed, sterile neutrinos can be produced in the Early Universe via active-sterile 
neutrino oscillations. Depending on the structure of the neutrino mixing matrix and the 
neutrino masses involved, sterile neutrino production can be copious. In this case, sterile 
neutrinos are expected to leave a clear imprint on cosmological observables, as active 
neutrinos do. Given the impressive accuracy reached by recent cosmological probes on 
one side, and the well-defined nature of the cosmological concordance model capable 
of interpreting those data on the other side, it is now possible to use the Universe to 
constrain sterile neutrino properties that have traditionally been probed via laboratory- 
based neutrino sources. 

As already extensively studied in the literature, we take into account how both 
neutrino interactions with the primeval medium, and neutrino coherence breaking 
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effects, affect the evolution of sterile neutrino abundances [121 EH EB] • Most previous 
work generally considers only the simplest sterile neutrino mixing matrix structures 
possible, that involves mixing a single sterile neutrino species with one or more active 
species |23j. In this paper we take a more general approach. We generalize those 
results to account for all mixing effects arising when multiple sterile neutrino flavors 
participate in the oscillations. Indeed, once sterile neutrinos are added to the Standard 
Model fermion content, there is no fundamental reason to expect only sterile single-flavor 
mixing effects. 

In addition to full numerical results that are valid, to a good approximation, for 
generic sterile neutrino models (with the only requirement of relativistic decoupling of 
sterile neutrinos), we try to emphasize the underlying physics by developing an analytical 
description that is applicable for sterile neutrinos that are significantly more massive 
than the active ones, for a mixing matrix that is approximately block-diagonal in its 
active and sterile sectors, and for sub-thermal sterile neutrino abundances at decoupling. 
In the latter case, we obtain the following three main results, valid for any number of 
sterile neutrino species: first, that the production of sterile neutrinos arises entirely from 
the breaking of the coherence of the evolution of active neutrinos; second, that each 
active flavor contributes independently to the sterile neutrino abundances; and third, 
that the evolution of the abundances of distinct sterile neutrino species is independent 
from one another, in absence of mixing in the sterile sector. 

Within our full multi-flavor framework, we have also studied in detail the 
phenomenology of the so-called (3 + 2) sterile neutrino models [301 El] • Two additional 
heavy (mostly-sterile) neutrino states with masses in the eV range, beyond the minimal 
three active neutrino mixing scenario, and with small amounts of electron and muon 
flavor content, are considered in this case. These (3 + 2) models were originally 
introduced as a means to reconcile via standard — ► v e oscillations the currently 
unexplained v e excess observed by the LSND experiment with the solar and atmospheric 
oscillation signatures. We first study the sterile neutrino evolution in the Early Universe 
of (3 + 2) models that offer a potentially viable explanation of all short-baseline 
oscillation data, including the LSND excess and first —>■ v e results recently released 
by the MiniBooNE Collaboration. We then contrast the expected energy density 
of the Universe in the form of sterile neutrinos obtained by those models with the 
latest cosmological observations coming from cosmic microwave background datasets 
and galaxy surveys, mostly (but not only) in the framework of the minimal ACDM 
cosmological model. 

In contrast with some literature, we find that fully solving the neutrino kinetic 
equations is necessary even for the relatively large active-sterile mixing implied by this 
class of models, since sterile neutrino states do not always feature thermal abundances 
at decoupling. This result, valid for the minimal ACDM cosmological model, should 
be taken into account in future studies on sterile neutrino production in the Early 
Universe, in order to avoid possibly misleading conclusions. We find that the (3 + 2) 
model that best describes short-baseline oscillation data is excluded at high confidence 
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level by cosmological observations, which do not allow for large amounts of energy 
density in the form of relativistic species in addition to the energy density of active 
neutrinos, unless non-minimal ACDM models are used. Nevertheless, by fully exploring 
the neutrino parameter space, we do find (3 + 2) sterile neutrino models that provide a 
perfectly acceptable description of short-baseline and cosmological data simultaneously, 
with a 11% probability for compatibility between the two datasets. As a consequence, 
we conclude that (3 + 2) models are significantly more disfavored by the internal 
inconsistencies between sterile neutrino interpretations of appearance and disappearance 
short-baseline data themselves [33] , rather than by the used cosmological data. Finally, 
a global analysis of all short-baseline plus cosmological data allow us to further constrain 
the (3 + 2) neutrino parameter space, in a region that should be accessible by future 
experiments aiming to study neutrinos, both the ones permeating the cosmos and 
neutrinos produced on Earth. 
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